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Abstract 



We compute the complete 0(a 2 s ) anomalous dimension matrix relevant for the b — > sj, 
b — >• sg and b —* s£ + £~ transitions in the standard model and some of its extensions. For 
radiative decays we confirm the results of Misiak and Miinz, and of Chetyrkin, Misiak 
and Miinz. The 0{a 2 s ) mixing of four-quark into semileptonic operators is instead a new 
result and represents one of the last missing ingredients of the next-to-next-to-leading-order 
analysis of rare semileptonic B meson decays. 



1 Introduction 



Inclusive radiative B decays place very important constraints on models of physics be- 
yond the Standard Model (SM). The present experimental world average for the branching 
fraction of the radiative B — > X s j decay is pQ 

BR(B -> X s7 ) cxp = (3.34 ± 0.38) x 10" 4 , (1) 

while the most recent SM prediction is [21131 

BR(B -> X s7 ) th = (3.70 ± 0.30) x 10~ 4 . (2) 

The experimental error is rapidly approaching the level of accuracy of the theoretical 
prediction. The main limiting factor on the theoretical side resides in the perturbative 
QCD calculation and is related to the ambiguity in the definition of the charm quark mass 
in some two-loop diagrams containing the charm quark To improve significantly the 
present Next-to-Leading-Order (NLO) QCD calculation, one would need to include one 
more order in the strong coupling expansion, and compute at least the dominant NNLO 
effects: a very challenging enterprise, which seems to have already captured the imagination 
of some theorists [3]. 

The present calculation of the branching ratio of B — > X s 7 consists of several parts that 
are worth recalling. Perturbative QCD effects play an important role, due to the presence of 
large logarithms L = lnrnj/M^, that can be resummed using the formalism of the operator 
product expansion and renormalization group techniques. The main components of the 
NLO calculation, which aims at resumming all the next-to- leading 0{a™L n ~ l ) logarithms, 
have been established more than six years ago. They are i) the 0(a s ) corrections to the 
relevant Wilson coefficients [El El El, H) the 0(a s ) matrix elements of the corresponding 
dimension-five and six operators |Hj, and Hi) the 0(a 2 s ) Anomalous Dimension Matrix 
(ADM) describing the mixing of physical dimension-five and six operators |HjllO | ITT | ll2j. 
After the 0(a s ) matrix elements of some suppressed operators have been computed last 
year [5J, the NLO calculation is now formally complete. Higher order electroweak [T3*1IT1] 
and non-perturbative effects ^ElEH] amount to a few percent in the total rate and seem 
to be well under control. However, as the actual measurements adopt a lower cut on the 
photon energy — typically around 2 GeV — the total rate is not accessible experimentally. 
In fact, unless the photon energy cut lies well below 2 GeV, the measured rate depends 
significantly on the non-perturbative structure of the B meson, and in particular on the 
Fermi motion of the bottom quark inside the B meson. The ensuing theoretical uncertainty 
is significant ^3] and is generally included in the reported experimental error, since the 
experimental measurements are extrapolated to conventionally defined total rates. 

Nearly all the ingredients of the NLO QCD calculation involve a considerable degree 
of technical sophistication and have been performed independently by at least two groups, 
sometimes using different methods. However, the most complex part of the whole enter- 
prise, the calculation of the two-loop dimension-five ^1] and of the three-loop dimension- 
six ^2] 0(a 2 s ) ADM, has never been checked by a different group. The main purpose of 
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this work is to present an independent calculation of the 0{a 2 ) ADM governing the b — > S7 
and b — > sg transitions. 

The rare semileptonic decay B — > X s i + i~ represents, for new physics searches, a route 
complementary to the radiative ones. Experimentally, the exclusive modes have been 
recently observed for the first time [IZj, and we now also have a measurement of the 
inclusive branching fraction ^Hj- More progress is expected in the future from the B 
factories. Because of the presence of large logarithms already at zeroth order in a s , a 
precise calculation of the B — * X s £ + £~ rate involves the resummation of formally next- 
to-next-to-leading 0(a™L n ~ 2 ) logarithms. The Next-to-Next-to-Leading-Order (NNLO) 
QCD calculation of B — > X s £ + £~ has required the computation of i) the 0(a s ) corrections 
to the corresponding Wilson coefficients [Zj and ii) the associated matrix elements at 
0(a s ) ^nH201- Moreover, it involves Hi) the 0{a 2 s ) ADM, but the operator basis must be 
enlarged to include the semileptonic operators characteristic of the b — > si + l~ transition. 
The 0(a> 2 ) mixing between four-quark and the relevant semileptonic operators is one of 
the last missing pieces of a full NNLO analysis of rare semileptonic B meson decays. Since 
we include the semileptonic operators in the 0(a 2 ) calculation, we are able to close this 
gap. Notice however, that a formally complete NNLO calculation of the B — > X s £ + £~ rate 
would also require the knowledge of the 0(a 2 ) self-mixing of the four-quark operators [21], 
of the 0(a s ) matrix elements of the QCD penguin operators, and of the renormalization- 
group-invariant two-loop matrix element of the vector-like semileptonic operator. 

We perform the calculation off-shell in an arbitrary gauge which allows us to ex- 
plicitly check the gauge-parameter independence of the mixing among physical operators. 
To distinguish between infrared (IR) and ultraviolet (UV) divergences we follow [22] arid 
introduce a common mass M for all fields, expanding all loop integrals in inverse powers 
of M. This makes the calculation of the UV divergences possible even at three loops, as 
M becomes the only relevant internal scale and three-loop tadpole integrals with a single 
non-zero mass are known. On the other hand, this procedure requires to take into account 
insertions of non-physical operators, as well as of appropriate counterterms. 

We have so far emphasized the inclusive modes, as they are amenable to a cleaner 
theoretical description. However, one should not underestimate the importance of the 
exclusive B meson decays like B -> K*^y [23, B -> K*£ + £~ [17]. B -» py [231I23| and B -> 
p£ + £~ . A thorough study of the exclusive channels can yield useful additional information in 
testing the flavor sector of the SM. These processes have received a lot of theoretical interest 
in recent years and their accurate description will also benefit from a firm understanding 
of higher order perturbative corrections. 

The ADM we have computed can be used in analyses of new physics models, provided 
they do not introduce new operators with respect to the SM. This applies, for example, to 
the case of two Higgs doublet models [21111112131212] , and to some supersymmetric scenarios 
with minimal flavor violation. See for instance [23121] • On the other hand, in left-right- 
symmetric models [2H1I2E! and in the general minimal supersymmetric SM [2H], additional 
operators with different chirality structures arise. In many cases one can exploit the chiral 
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invariance of QCD and use the same ADM, but in general an extended basis is required. 

Our paper is organized as follows: in Section 2 we recall the relevant effective Lagrangian 
and list all the dimension-five and six operators that will be needed in the calculation. In 
Section 3 we review the renormalization procedure and explain how the operator renor- 
malization matrix can be extracted from the matrix elements at higher orders. The actual 
two- and three-loop calculation is described in Section 4, while the results for the ADM 
are presented in Section 5. We provide some intermediate results which can be useful in a 
number of related applications in an Appendix. They include the complete 0(a s ) and the 
relevant entries of the 0{a 2 s ) operator renormalization matrices. 

2 The Effective Lagrangian 

Let us briefly recall the formalism. We work in the framework of an effective low-energy 
theory with five active quarks, three active leptons, photons and gluons, obtained by 
integrating out heavy degrees of freedom characterized by a mass scale M > M w . In the 
leading order of the operator product expansion the effective off-shell Lagrangian relevant 
for the b — > S7, b — > sg and b — > s£ + £~ transition at a scale \x is given by 



Here the first term is the conventional QCD-QED Lagrangian for the light SM particles. 
In the second term Vij denotes the elements of the CKM matrix and Cj(/x) are the Wilson 
coefficients of the corresponding operators Qi built out of the light fields. 

In our case it is useful to divide the local operators Qi entering the effective Lagrangian 
into five different classes: i) physical operators, ii) gauge-invariant operators that vanish by 
use of the QCDxQED Equations Of Motion (EOM), Hi) gauge- variant EOM- vanishing op- 
erators, and iv) so-called evanescent operators that vanish algebraically in four dimensions. 
In principle, one could also encounter v) non-physical counterterms that can be written as 
a Becchi-Rouet-Stora-Tyutin (BRST) variation of some other operators, so-called BRST- 
exact operators. However, they turn out to be unnecessary in the case of the 0{a 2 s ) mixing 
of the operators Qi considered below. See also jTTj. 

Neglecting the mass of the strange quark and the electroweak penguin operators which 
first arise at 0(a), the physical operators can be written as |3U] 



C-cfi — £qcdxqed(^, d, 



s, c, b, e, ji, t) + — 7=-V t * s V tb C iW Qi ■ 



(3) 



i=X 



Qi = (s Ll ,T a c L )(c L ^T a b L ), 
Q2 = (sl7/.cl)(c L 7 m 6 l ) , 

Qs = (sLlnb L ) V (qYq) , 




Q5 = {sLl^lulpb L ) V [qi^l u l p q) , 
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-^m b (s L a tlv b R )F IJV , 
9 2 

Qs = -m b (s L a^T a b R )Gl 



Q 



Q 9 = ^(^7m^)V {hn), 

9 1 

' s L i»b L )Y,Xh^ 
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g 



(4) 



where the sum over q and I extends over all light quark and lepton fields, respectively, 
e (g) is the electromagnetic (strong) coupling constant, q L and q R are the chiral quark 
fields, F^y (C" ) is the electromagnetic (gluonic) field strength tensor, and T a are the color 
matrices, normalized so that Tr (T a T b ) = S ab /2. Notice that since QCD is flavor blind it 
is not necessary for our purposes to consider the analogues of Q\ and Q2 involving the up 
instead of the charm quark. 

The physical operators given above consist of the current-current operators Qi and 
Q 2 , the QCD penguin operators Q3-Q6, the magnetic moment type operators Q 7 and Q 8 , 
and the semileptonic operators Qg and Q10, relevant for the b —> s£ + £~ transition. We 
have defined Qi-Qe in such a way that problems connected with the treatment of 75 in 
n = 4 — 2e dimensions do not arise |10j . Consequently, we are allowed to consistently use 
fully anticommuting 75 in dimensional regularization throughout the calculation. 

The gauge-invariant EOM-vanishing operators can be chosen to be jHllZj 



Qn 
Q12 
Q13 

Ql4 

Q15 
Qw 



- 2 s Ll ^F, v + - 2 {s L l,b L ) V Q/t/YV) 



9 

-s L ^T a b L D»G a uu + Q 4 



9 



9' 

% 

9 

ie 

IF 

% 

9 



m h s L ppb R , 
2 s L pppb L , 

s L p a^b L F^ - F^s^pbL 



(5) 



+ Q 



7 ■ 



s L p a^T a b L Gl v - G a -s L T a a^pb L 



where the sum over / runs over all light fermion fields, while and denotes the 
covariant derivative of the gauge group SU(3)c x U{1)q acting on the fields to the right 
and left, respectively. Notice that the set of operators Qi-Q\q closes off-shell under QCD 
renormalization, up to evanescent operators |30|l7j. 
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In order to remove the divergences of all possible one-particle irreducible (1PI) Green's 
functions with single insertion of Qi—Qio we also have to introduce the following gauge- 
variant EOM-vanishing operators 



Q 17 = -m b s L 
9 



pG-Gp 



>R 



Ql8 = i 



s L (pGG- GGp b L - im b s L GGb R 



Q 



Q20 = i 



Q21 
Q22 
Q23 



s L [pp G + Gppj b L + im b s L @pb R 
s L (p G®G afl - G^G^pj b L - im b s L GlG a »b R 
s L (pD^ G» + G^pj b L + im b s L G^b R 



(6) 
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Q24 = d 



abc 



s L lp T a + T a pj b L + im b s L T a b R ^ d^Gl , 
s L p Gpb L + im b s L p Gb R , 
s L (p T a - T a p] b L - im b s L T a b R 



G\G C » , 



where G° denotes the gluon field, and we have used the abbreviations G^ = G®T a and 
d abc = 2Tr({T a , T b }T c ). 

It is important to remark that the EOM-vanishing operators introduced in Eqs. (JHJ) and 
(JHJ) arise as counterterms independently of what kind of IR regularization is adopted in the 
computation. However, if the regularization respects the underlying symmetry, and all the 
diagrams are calculated without expansion in the external momenta, non-physical operators 
have vanishing matrix elements [3T]. In this case the EOM-vanishing operators given in 
Eqs. © and © play no role in the calculation of the mixing of physical operators. If the 
gauge symmetry is broken this is no longer the diagrams with insertions of non- 

physical operators will generally have non-vanishing projection on the physical operators. 
Since our IR regularization implies a massive gluon propagator, non-physical counterterms 
play a crucial role at intermediate stages of the calculation. 

As far as the evanescent operators are concerned, another eight operators are needed 
in order to find the 0{a 2 s ) mixing of the physical operators Q\-Qiq. Following [T21ITU] they 
can be defined as 

Q25 = (s L l,lul P T a c L )(c L YYl p T a b L ) - 16Q 1 , 

Q26 = {s L ^^ u - 1p c L ){c L Yl u l p b L ) - I6Q2 , 

Q27 = {sLl^lplvlM Y] (qYYYYYq) + 64Q 3 - 20Q 5 , 
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Q28 = (s Ll ^ ulplalr T a b L ) ^ 9 ( g y 7 Y7Yr a g) + 64Q 4 - 2oq 6 , 

Q29 = (s Ll ^ lplalT T a c L )(c L rrYl a YT a b L ) ~ 256Q! - 20Q 25 , 

Q30 = (sL7/,7,7p7 CT 7rCL)(c L7 At 7 l '7 P 7 CT 7^L) - 256Q 2 - 20Q 26 , 

Qsi = {sLl^ulpl.lrlvlM ^JgYT^YTYT^) + 1280Q 3 - 336Q 5 , 

Q32 = (sL7/,7,7p7a7r7^ a 6 L ) J2 (q^YYYYYl^q) + 1280Q 4 - 336Q 6 . (7) 

3 Renormalization of the Effective Theory 

Our aim is to study the renormalization properties of the physical operators Q1-Q10 in- 
troduced in Eq. (J3J). Upon renormalization, the bare Wilson coefficients C^sip) of Eq. (jSJ) 
transform as 

C i , B {ii)=Z ji Cj{ i i), (8) 
where the renormalization constants can be expanded in powers of a s as 

fe=l ;=o 

Following the standard MS scheme prescription, is given by pure 1/e poles, except 
when z = 25-32 and j 7^ 25-32. In the latter case, the renormalization constant is finite, to 
make sure that the matrix elements of the evanescent operators vanish in four dimensions 
|32U33j . The calculation of an effective amplitude A c s, also involves the matrix element 
(Qi) = (F\Qi(fi)\I) of the operator Qi between a initial state / and a final state F, 
which is renormalized by the usual coupling, mass and wave function renormalization 
factor characteristic of the operator, Qi — > Z(Q{). The renormalized effective amplitude is 
therefore given by 

A fS = Z ji C j { i i)(Z{Q i )) R , (10) 

where (Z(Qi)) R denotes the matrix element of the operator Z(Qi) after performing cou- 
pling, mass and wave function renormalization. Clearly, it is also possible to define the 
operator renormalization constant Z%j from the relation between unrenormalized and ampu- 
tated Green's functions via (Z(Qi)) R =Zij{Qj) B . In this case, one simply has-Z^- = Z^ . 
In general Z(Qi) will not be proportional to Qi. For example, in many of the EOM- 
vanishing operator introduced in Eqs. (J^J) and (jUJ) one has two different terms, only one of 
which has a factor of m;,. Correspondingly, the mj renormalization of the operator is 

Z mb (Q l )=Q i + (Z mb -l)Q' l , (11) 

where Z mb denotes the mass renormalization constant of the bottom quark, and Q\ is the 
part of Qi proportional to m^. 
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Figure 1: Some of the two- loop 1PI diagrams we had to calculate in order to find the 
0{a 2 g) mixing of the complete set of operators Qi-Qm- 



The product on the right-hand side of Eq. (jlOJ) must be finite by definition at any given 

•(*: 

i.i 



order in a s . Therefore, requiring the cancellation of UV divergences we can extract Z,' 



order by order. The result, up to third order, reads 

zS ) (Qj) ( k ) --{z(Qi)) ( k\ 

z$(Qi)f = -<z(Qi))$ - z% ] (z(Qi))R , (12) 

ZifiQj^R = -(Z(Qi))f - z\f{Z{Q, j ))^ - zgiZiQi))® , 

where the superscript (k) always stands for the k-th order contribution in a s . 

If we leave aside the complication that in general Z(Qi) will not be proportional to Qi, 
and write symbolically (Z(Qi))n = Z^Q^b the above relations can be rewritten in terms 
of bare quantities. Up to third order in a s we obtain 

Zij\Qj)s = ~(Qi) < h > ~ Zi(Qi) B , 

- Z$ zf\ Qj )f - Z?\Q,)f , 

- zfzf m f - z$zf\ Qi )f - zfmf ■ 

The first line in Eqs. f)13j) recalls the familiar result that the one-loop renormalization 
matrix is given by the UV divergences of the one-loop matrix elements, after performing 
wave function and possibly coupling and mass renormalization. For example, in the case of 
the operators Qi-Qe, one has Zj = Z^ with Z q denoting the wave function renormalization 
constant of the quark fields, and Eqs. (|TH|) take a particularly simple form, which upon 
expansion in a s reproduces the classical results derived more than ten years ago [52] . 

For a given set of operators and knowing the QCD renormalization constants, the solu- 
tion of the above systems of linear equations requires the calculation of a sufficient number 
of Green's functions for different external fields with single insertions of the operators Qi. 
In our case, in order to determine the complete Z^ of all the operators introduced in 



(13) 



S 



Figure 2: Some of the three-loop 1PI diagrams we had to calculate in order to find the 
mixing of the four-quark operators into Q7-Q10 at 0(af). 



Section 2, it is sufficient to calculate the 0(a k s ) matrix elements of Q1-Q32 for the b — ■> sec, 
b — > s, b — > sj, b — > sg and b — > sgg transitions — see Fig. ^ As we are interested in a 
subset of the three-loop ADM, we have actually calculated only the three-loop b — > S7 and 
b ^ sg amplitudes involving insertions of — see Fig. |3 We have calculated the 

complete off-shell amplitudes up to terms proportional to external momenta squared. By 
using the EOM it is therefore straightforward to extract the mixing into Q7-Q1Q. Notice 

(k) 

that the results for the cannot depend on the considered Green's functions and that the 
pole parts need to have the structure of the complete set of local operators Qi-Qw Both 
features represent a powerful consistency check of the computation of the renormalization 
constants . 

The normalization of the physical operators adopted in Section 2 has been chosen 
[7j in such a way that the power of a s in Z^ is equal to the number of loops of the 
contributing diagrams. For instance, without the factor 1/g 2 in Q7-Q10, as in the standard 
normalization adopted in PUSlEOj, both one- and two-loop diagrams contribute to the 
0(a s ) mixing matrix, because of the 0(a s ) two-loop mixing of four-quark into magnetic 
operators. This choice simplifies both the implementation of the renormalization program 
and the resummation of large logarithms, since the redefinition enables one to proceed for 
b — > s£ + £~ in the same way as in the 6 — » 57 and b — > sg case. 

In a mass independent renormalization scheme Z^ is /i-independent. This allows to 
check the renormalization of two- and three-loop matrix elements. The right-hand sides of 
the Eqs. (112)) and (|13|) receive contributions from irreducible two- and three-loop diagrams 
as well as one- and two-loop counterterms. The /i-dependence is different in each case and 
governed by the n-loop factor (fi 2e ) n . The UV structure of the k-th term is therefore given 
by 

k n 1 

^^^EEW"/^. (14) 

n=0 1=1 

where M^ 1 ^ denotes the l/e l pole of the sum of all n-loop contributions. Expanding 
in powers of e we find the following set of equations which have to be fulfilled to get a 
//-independent Z^ up to three-loop order: 

3M^ +2M^ + =0, 
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3M (3 ' 3) + 2M (2 ' 3) + M (1 - 3) = , 
9M (3 ' 3) + 4M (2 ' 3) + M (1 - 3) = . 



(15) 



This system of equations provides us with a powerful check of the renormalization of two- 
as well as three-loop diagrams. Notice that the locality of UV divergences also places some 
constraints on the renormalization matrix itself. We will turn to this point later on. 

4 The Calculation of the Operator Mixing 

In the renormalization of QCD and QED at higher orders the standard method of extracting 
the UV divergence structure of a Feynman integral is to perform the calculation with 
massless propagators. However, if one uses massless propagators to compute three-point 
or higher Green's functions one might generate spurious IR infinities which, in dimensional 
regularization, cannot be distinguished from the UV divergences one seeks. There exist 
several methods [3^ to overcome this problem, but they are generally quite involved and 
not suitable to the automated evaluation of a large number of diagrams. In the approach 
of [22] the so-called IR rearrangement is performed by introducing an artificial mass. For 
the calculation of the renormalization constants this means that we can safely apply Taylor 
expansion in the external momenta after introducing a non-zero auxiliary mass M for each 
internal propagator, including those of the massless vector particles. The auxiliary mass 
regulates all IR divergences and the renormalization constants can be extracted from the 
UV divergences of massive, one-scale tadpole diagrams, that are known up to the four-loop 
level 

Following references f2E2] , the starting point of our procedure is the exact decompo- 
sition of a propagator: 

1 _ 1 p 2 + 2k ■ p - m 2 + M 2 1 

(k + pf -m?~ k 2 -M 2 k 2 -M 2 (jfc + pf - m * ' 

Here k is a linear combination of the integration momenta, p stands for a linear combination 
of the external momenta, and m denotes the mass of the propagating particle. If we assume 
that the dimensionality of the operators in our effective theory is bounded from above, and 
we apply recursively the above decomposition a sufficient number of times, we will reach the 
point where the overall degree of divergence of a certain diagram would become negative 
if any of its propagators were replaced by the last term in the decomposition. We are then 
allowed to drop the last term in the propagator decomposition, as it does not affect the 
UV divergent part of the Green's function after subtraction of all sub divergences. 

This algorithm can be also simplified by the following observation [22] • The terms 
containing powers of the auxiliary mass squared in the numerators contribute only to UV 
divergences that are proportional to those powers of M 2 . The latter are local after the 
subtraction of all subdivergences, and must precisely cancel similar terms originating from 
integrals with no auxiliary mass in the numerators. Since the decomposition of Eq. ()16|) 
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is exact, no dependence on M 2 can remain after performing the whole calculation. This 
observation allows one to avoid calculating integrals that contain an artificial mass in the 
numerator. Instead of calculating them, one can just replace them by local counterterms 
proportional to M 2 which cancel the corresponding sub divergences in the integrals with 
no M 2 in the propagator numerators. Nevertheless, the final result for the UV divergent 
parts of the Green's functions are precisely the same as if the full propagators were used. 

The counterterms proportional to M 2 in general do not preserve the symmetry of the 
underlying theory, specifically they do not have to be gauge-invariant. Fortunately, the 
number of these counterterms is usually rather small, because their dimension must be two 
units less than the maximal dimension of the operators belonging to the effective theory. 
For instance, in QCD only a single possible gauge-variant operator exists that fulfills the 
above requirement. It looks like a gluon mass counterterm, 



and cancels gauge-variant pieces of integrals with no M 2 in the numerators. To ensure 
that our renormalization procedure with the fictitious gluon and photon mass is valid, we 
have checked explicitly the full MS renormalization of QCD and QED up to the three-loop 
level, finding perfect agreement with the results given in the literature jHEj- In our case, 
beside the term in Eq. (117)1 . we also have M 2 counterterms of dimension-three and four, 
some of which explicitly break gauge invariance: 



where denotes the photon field. 

As already mentioned in Section 2, another side effect of our IR regularization is that 
we have to consider insertions of non-physical effective operators in our calculation. Let 
us explain this point in more detail. Non-physical counterterms generally arise in QCD 
calculations, but the projections of their matrix elements on physical operators vanish un- 
less the underlying symmetry is broken at some stage. Due to the exact nature of the 
decomposition Eq. ()16|) . the UV poles of the diagrams obtained by our method are correct 
after the subtraction of all sub divergences. However, the UV poles related to subdiver- 
gences and their subtraction terms both depend on the finite parts of certain lower loop 
diagrams, which in our approach are not necessarily correct and do not comply with the 
usual Slavnov- Taylor identities. For instance, the introduction of the IR regulator invali- 
dates the argument that guarantees vanishing on-shell matrix elements for the non-physical 
operators. One therefore expects non-negligible contributions to the counterterms from all 
possible operators with appropriate dimension. Consequently, all EOM-vanishing opera- 
tors, gauge- invariant or not, and in general even BRST-exact operators must be included 
in the operator basis. The "incorrect" sub divergences are present in both counterterm and 
irreducible diagrams, but they cancel in their sum, provided the calculation is carried out 
in exactly the same way. The operator renormalization constants calculated in this way 
are correct for all the operators in the complete basis. 



M 2 G a ..G a ^ , 



(17) 



g 2 



M 2 
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The large number of diagrams which occurs at higher orders makes it necessary to 
generate the diagrams automatically. For the evaluation of the ADM presented here all 
diagrams have been generated by the Mathematica jHZj package FeynArts jHHj, which 
provides the possibility to implement the Feynman rules for different Lagrangians in a 
simple way. We have adapted it to include the effective vertices induced by the operators 
Qi-Qzi- We have processed the FeynArts output using two independent programs. In 
one case the output is converted into a format recognizable by the language Form [3*§] . 
The group theory for each graph as well as the projection onto all possible form factors 
is performed before the integrals are evaluated. The very computation of the integrals is 
done with the program package MATAD |3Q] , which is able to deal with vacuum diagrams at 
one-, two- and three-loop level where several of the internal lines may have a common mass. 
The calculation of the tadpole integrals in MATAD is based on the so-called integration-by- 
parts technique jlj. The second program is entirely a Mathematica code, which for the 
three-loop integrals uses the algorithm described in detail in [2*2] . 



5 The Anomalous Dimension Matrix 



The anomalous dimensions 7^ defined by 



can be expressed in terms of the entries of the renormalization matrix Zjj as follows 

d_ 



7y = Zikfi—Z^ 1 . (20) 



In a mass independent renormalization scheme the only /i-dependence of resides in the 
coupling constant. In consequence, we might rewrite Eq. ()20|) as 

ll3 = 2P(e,a s )Z tk ^-Z~l, (21) 

where /3(e, a s ) is related to the f3 function via 

(3{e,a s ) = a s (-e + {3{a s )) . (22) 

The finite parts of Eq. (|2*T|) in the limit of e going to zero give the anomalous dimensions. 
Expanding the anomalous dimensions and the (3 function in powers of a s as 



00 , 00 , 

and «"•> = -£(£) (23) 

k=l k=l 
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we find in accordance with 



up to third order in a s : 



T 
7 
1 



(1) = 4^(2,i) _ 2Z<MzW> - 2Z^Z^ + 2/3b# 1,0) , 

( 2 ) = 6Z (3>1) - 4Z (2,1) Z (1 ' 0) - 2Z (1 ' 1} Z {2 ' 0) - 4ZP' 0) Z (1 '^ - 2Z (1 ' 0) Z (2,1) 
+ 2Z^Z^Z^ + 2Z^Z^Z^ + 2Z^Z^Z^ 
+ 2/3 1 Z( 1 -°) + 4/? Z( 2 '°)-2/3 Z( 1 ' ^( 1 ' ) . 



(24) 



On the other hand the pole parts of Eq. (|21J) must vanish. From this condition one obtains 
relations between single, double and triple 1/e poles of the Zy, which constitute a useful 
check of the calculation. In agreement with [22] we find 

2(2,2) = 1^(1,1)^(1,1) _ I^d.l) ; 



^(3,3) = l _z^Z^Z^ - -foZ^zW + , 

2(3,2) = 22(2,i)2(i,i) i i2 (1 ' 1) z (2,1) - i2 (1 ' 1) z (1 ' 0) z (1,1) - iz (1,0) i' (1 ' 1) z (1 ' 1) 

3 3 3 6 

3 3 o 



(25) 



Having summarized the general formalism and our method, we will now present our 
results for five active quark flavors. For completeness we start with the regularization- and 
renormalization-scheme independent matrix which is given by 



1 



(o) 



/ "4 f 
12 










\ 













4 

o -m o 



ill 

!) 



2'.« 

9 








2 
" 9 
4 
3 

HI 

:! 
11)1) 
9 

:', 

!) 











iii 

!) 










2 

5_ 

6 

20 

_ 2 
3 

















14 
:s 

!) 
















-6 





_32 
27 

_8 

!) 

_!<! 

9 

'HI 

27 
112 
9 

512 
27 








ill 













_46 
3 



\ 



/ 



(26) 



While the matrix ■y^ is renormalization-scheme independent, 7 ( - 1 ^ and are not. In 
the MS scheme supplemented by the definition of evanescent operators given in Eq. (J7J) we 
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obtain 



/ _355 _ 502 1412 1369 134 35 232 167 2272 

9 97 243 243 243 162 243 162 729 



9 
_ 35 
3 







~.l>_i 

27 
28 
3 









1412 
243 
416 

81 
4468 

81 
8158 
243 
25168Q 

81 



1369 
243 
1280 
81 

31469 

81 
59399 

243 
128648 



134 
243 
56 
Si 
4110 
81 
269 
486 
23836 
81 



35 
162 
35 
27 
3373 
108 
12899 
648 
6106 
27 



232 
" 243 
464 
81 

<±k 
81 

200 
243 
6464 
81 



167 
162 
76 
27 
36.N 
27 
1409 
162 
13052 
27 



2272 
729 
1952 
243 
6752 
243 
2192 
729 



and 



7 



(2) 













58640 
243 




26348 
243 


14324 
243 


2551 
162 


11408 2740 37856 
243 81 729 


























2600 


























2192 1975 f) 
81 27 
























































/ 


? 


? 


? 


? 


? 


? 


13234 
2187 


13957 
2916 


1359190 | 6976 a 
19683 ~ l ~ 243 ^ 3 




? 


? 


? 


? 


? 


? 


20204 
729 


14881 
972 


229696 3584 a 
6561 81 * 3 










? 


? 


? 


? 


92224 
729 


66068 
243 


1290092 , 3200 a 
6561 ^ 81 1,3 










? 


? 


? 


? _ 


184190 
2187 


1417901 
5832 


819971 19936 a 
19683 243 ^ 3 










? 


? 


? 


? 


1571264 


3076372 


16821944 | 30464 a 












729 


243 


6561 1 81 *> 3 










? 


? 


? 


? _ 


1792768 
2187 


3029846 
729 


17787368 286720 a 
19683 243 ^ 3 






















? 




























? 


? 































9769 
27 


\ 






























o \ 












232 
3 



(27) 



7 



o \ 

o 











27 



(2? 



/ 



The question marks stand for unknown entries that we have not calculated. As far as the 
remaining entries are concerned, our results for the one- and two-loop mixing of Q\-Q§ 
agree with those of ^H], and therefore also with previous results jHj that were obtained 
in a different operator basis [12! • We also confirm other well-established results, like the 
two- loop mixing of Qi~Qe into Q7-Q10 [12] • See [7] for the conversion to our basis. For 

the two- loop mixing of Q7 and Q 8 we confirm the results of [TT] . and for the entries 7^ 

(2) 

and 7§j , that is the three-loop mixing of four-quark operators into magnetic ones, those 
of ^2]. Neither of the latter results have been checked before, with the exception of the 
self- mixing of Q7 [HJ. 

We recall that the ADM reported in ^2] refers to so-called effective coefficients. The 

(2) 

relation with our results will be discussed below. Finally, the entries 7^ which contain 
terms proportional to the Riemann zeta function £3, are entirely new. The ADM entries 
involving non-physical operators can be readily obtained from the renormalization matrices 
provided in the Appendix. Finally, let us mention here once again that for a NNLO analysis 
of rare semileptonic B decays one would also need to know the self-mixing of Qi~Qe, which 
will be discussed in a separate communication [21]. 

It might be useful to recall explicitly the relation between the ADM in our basis and the 
ADM in an operator basis where Q7-Q10 are not rescaled by 1/g 2 . The latter is frequently 
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used for phenomenological applications 
Ci(fi), are given by 

CAu) , for % 



The Wilson coefficients in that basis, 



1-6, 
7-10 



(29) 



while the coefficients in the expansion in powers of a s of the corresponding anomalous 
dimensions 7^, take the following form: 



~(fe-i) 

lij 



(fe-1) 
lij 

(k) 
lij 5 



y ■ 



for i,j = 1-6 , 
for i = 1-6 , and j 
for i,j = 7-10 . 



7-10, 



(30) 



Notice that the perturbative expansion of the anomalous dimensions 7^, that govern the 
evolution of Ci(n) already starts at zeroth order in a s , whereas 7^ contains no such terms. 

In the case of radiative B decays it is sometimes convenient to use effective Wilson 
coefficients C 4 cff (/i) [IHj, defined in such a way that the leading order b — > 57 and b — > sg 
matrix elements are proportional to Cf (//) and C| ff (/i), respectively. In particular, drop- 
ping the semileptonic operators which are irrelevant for radiative B decays, the effective 
Wilson coefficients in the operator basis of ^2] are given by 



cm 



for % 



1-6, 
7-8. 



80^ 
9 ' 



(31) 



and 



In the MS scheme with fully anticommuting 75 one has = (0, 0, — |, — |, — y , — - 
y( 8 ) = (0, 0, 1, — |, 20, — y). In consequence, the first two coefficients in the perturbative 
expansion of the effective ADM 7 cfr , that governs the evolution of Cf (/x) read 

\ 



7' 



cff(0) 



/ 


-4 


8 
3 





2 
9 








208 
243 


173 
162 




12 








4 

3 








416 
81 


70 
127 













52 
3 





2 


176 
81 


14 
127 










40 


100 


4 


5 


152 


587 




9 


9 


!) 


(i 


243 


162 













256 
3 





20 


6272 
81 


6596 
27 










256 


r,(, 


40 


•2 


4624 


4772 




9 


9 


9 


■i 


243 


81 






















:',2 
3 

























32 
9 


28 
3 



(32) 



15 



and 



/ _35J> _ 502 1412 1369 134 35 818 3779 \ 

9 27 243 243 243 162 243 324 



7 



cff(l) 



355 
9 


502 
27 


1412 
243 


1369 
243 


134 
243 


35 
162 


818 
243 


377!) 
324 


35 
3 


28 
3 


416 
81 


1280 
81 


56 
81 


35 
27 


508 
81 


1841 
108 








4468 


31469 


400 


3373 


22348 


10178 


81 


81 


Si 


108 


243 


81 








8158 


59399 


20!) 


12899 


17584 


172471 


243 


243 


486 


648 


243 


648 








251680 
81 


128648 
81 


23836 
81 


6106 
27 


1183696 
729 


2901296 
243 








58640 
243 


26348 
243 


14324 
243 


2551 
162 


2480344 
2187 


3296257 
729 




















4688 
27 























2192 
81 


4063 
27 



(33) 



which are in perfect agreement with 



6 Summary 

We have recalculated the complete 0{a 2 s ) mixing relevant for the NLO analysis of radiative 
B decays in the SM and some of its extensions, confirming all previous results. In addition, 
we have also calculated the 0{a 2 s ) three-loop mixing between four-quark and semileptonic 
operators relevant for the NNLO analysis of rare semileptonic b — > s£ + £~ transitions. This 
is a new result and provides one of the last missing pieces in the NNLO calculation for 
B — > X s £ + £~ . We plan to complete the calculation of the NNLO ADM required for this 
process and to study its phenomenological implications in a separate communication. 

The calculation involves the UV divergences of diagrams up to three loops. Our results 
have been subject to several cross-checks, following from: i) the locality of the UV diver- 
gences, ii) the independence of the ADM from the external states used in the calculation, 
Hi) the completeness of our operator basis, iv ) the gauge-parameter independence of the 
mixing among physical operators, and v) the absence of mixing of non-physical into phys- 
ical operators. We have also reproduced the full MS renormalization of QCD and QED up 
to the three-loop level. 
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Appendix 



A.l One- and Two-Loop QCD Renormalization Constants 

The renormalization of the conventional QCD Lagrangian containing a massive bottom 
quark proceeds as usual. First, we introduce the renormalized fields and variables via 

G%B = 4 /2 ^ , u| = ^ V , q B = Zfq , 
g B = Z g g , m b)B = Z mb m b , M B = Z M M , 

where the subscript B denotes bare quantities and u a are the ghost fields. The gauge- 
parameter £ is kept unrenormalized. This is legitimate, because the non-renormalization of 
the gauge-parameter is guaranteed by the usual Slavnov- Taylor identity, which is unaffected 
by the IR regularization adopted for the Yang-Mills theory. On the other hand, our IR 
rearrangement requires the introduction of the gauge- variant subtraction in Eq. ([17jL which 
can be interpreted as a counterterm for a fictitious gluon mass M. 

Using the notation introduced in Eq. (JD*)). the MS renormalization constants at one-loop 
order take the following form 

z »' i)= (i-Kk 

Z ' A = -tc F , (35) 

= ~3C F , 

where Ca = 3 and Cp = 4/3 are the quadratic Casimir operators of SU(3)c- As usual 
Nf stands for the number of active quark flavors. Our result for Z^' 1 ^ agrees with the 
expression given in [22] . 

At the two-loop level the poles of the MS renormalization constants are given by 



7 (2,1) 


/59 
\16 


zt 1] = 


/95 
V96 




3^,2 

4° F 



25 . 1 i2 \ 



e ) C F C A + -C F N f 
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and 



Z^ = -^Cl + ^C F N f + °-C A N f , 
Z^--^ F -%C F C A + 5 -C F N f , 



^; 2) = jjc£ + yC f c a - c F Nf , 



(37) 



M V 384 192 s 128 V 2 V 12 24 7 ; 3 J 

Except for Z^' 1 ^ and Z^' , which have never been given explicitly, our renormalization 



constants agree with the results in the literature |4b], if one bears in mind that the original 
papers contain some typing errors. We have also calculated the three-loop renormalization 
constants [SH], but we do not report them here, as they are not needed in our calculation. 



A. 2 The Complete Operator Renormalization Matrix 

The general structure of the operator renormalization matrix is 



/ 



V 



A(fc,J) 
Zj pp 

A(k,l) 

n NP 



'EP 



A(k,J) 
PN 
A(k,l) 
^NN 



KM) Mk,l) 



'EN 



NE 



r (k,l) 
'EE 



(3« 



where P = 1-10 denotes the physical operators, N = 11-24 the EOM- vanishing operators, 
and E = 25-32 the evanescent operators. Throughout this section we set Nf = 5. 

The mixing of non-physical into physical operators must vanish at all orders in a s . This 
is in fact only a requirement on the ADM, but we have seen in Eq. (|24j) that the one-loop 
renormalization matrix Z^ 1 ' 1 ' is proportional to 7^, and therefore at one-loop it implies 
the vanishing of Z^'p^ and Z Fp ■ Since 7^) for the physical operators can be found in 
Eq. (|26|). it is sufficient to give here only the non-physical parts of Z 1 ^ 1 ^ and Z^ l,1 \ By 
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definition, the only non- vanishing parts of Z^ 1,0 ^ are Z^p^ and Z^ff- We find 



>(i.o) 



and 



64 
48 


32 
3 

-64 






4 
9 

3 


























04 
27 
Hi 
!) 


o\ 










8960 
3 


-2432 




1280 
3 


320 


64 
3 


-64 


16 







3840 



640 


4480 
9 




9464 
3 

16 




640 

!) 






1280 
3 




256 
9 




32 
3 






256 
3 

256 
■', 






2880 - 


3840 





-96 






















64 











609280 
3 


-16076 




98560 
3 


24640 


512 

3 


-512 


544 











304640 


630256 




49280 


98560 


2048 


256 




11264 fl I 


9 


3 




9 




3 




9 


3 




3 


" / 




/ 64 

27 


4 
9 





























o\ 






16 
9 


8 




































16 


192 


































(1,0) 


256 
3 


168 


































EN — 


256 
3 


-16 




































64 


96 




































544 


8448 




































i 11264 
\ 3 


6992 





























o J 





(39) 



(40) 



The 6x4 block in the upper left corner of Z^ p agrees with the expression for the upper 
6x4 block of c given in Eq. (46) of [TO] . 

The one-loop mixing of physical into non-physical operators is described by Z p 1 ^ and 
Z%£\ We get 



J PN 



( 


16 
27 


1 

!) 



































o\ 




4 
9 


2 
:) 








































8 
9 


4 
■ >, 








































16 
27 


28 
9 








































56 
9 


64 
3 








































256 
27 


268 
9 


























































































-8 











9 
4 



































































V 












































(41) 



19 



and 



KM) 

J PE 



/ 


5 

112 


2 

!) 


o 


o 


o 


o 


o 


\ 




1 






















































































1 






















2 

1) 


5 

12 




























































































{ 























o J 



(42) 



The 4x6 block in the upper left corner of Z p 1 ^ ag: 
block in the upper left corner of b given in Eq. (45) 
At one-loop we have moreover the mixing amon; 



rees with the expression for the 4x6 

of [ID]. 

EOM-vanishing operators, given by 



J NN 



( 


23 
3 






































o \ 







!) 
2 





4 

3 





■', 

H 


!) 
8 


9 
8 


3 
2 


3 

16 


3 
2 





!) 

4 


9 
16 







- 


11 4 c 

3 3« 











































4 


23 4 c 
3 3 ? 


























f+Ic 



















23 
3 





































-8 


8 


8 
9 


17 

4 





3 
4 





1 

8 








9 
2 


3 
8 










3 s 











11 

3 








































97 3 c 
24 8 ? 





31 
12 


13 
3 


61 It 
72 8 s 


25 
3 


-4 


17 17 c 

2 6 ? 


95, 3 c 
48 16 ? 










-4-fe 







49 , 3 c 
24~ r 8 s 





13 
12 


-8 


72 ' 8 5 


1 

3 





2^6 S 


48 16 s 



















1 

2 








1 


65 3 c 
12 4 s 


-1 


-3 


-3-£ 













-u 




2 
9 


17 

24 





1 

12 


7 

24 


144 ^8? 


13 

2 


1 

4 


4 12 s 


29 3 c 
48 16 s 













-£ 

3 s 

















"ft 





65 
12 


4 s 


16 s 








































65 7 c 
12 12 s 





{ 

















5 

12 





5 

12 


5 
(i 


5 

72 


10 

3 


5 
2 


5.5c 
2 T 6 S 


_«_3c j 
8 8 s / 



(43) 



and the mixing among evanescent operators, which reads 



/ 



'EE 



-7 


4 

3 





-6 














64 
3 








28 
9 


























1792 
3 








1568 
9 






-14 

13 
3 




-784 

2212 
3 

Ki.i) 



5 


2 


12 


<) 


1 

















13 


28 


3 


9 


-14 


64 
3 



















2 
!) 






o \ 



1 

12 






-64 38 

76 166 / 

9 3 / 



(44) 



The 4x4 block in the upper left corner of Z E ' E agrees with the expression for the 4x4 
block in the upper left corner of d given in Eq. (47) of |10j . The last block Z^ E \ contains 
only zeros. 
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Now we can proceed to the two-loop renormalization matrices. The non-vanishing 
blocks of Z^ 2 ' ^ are Z^p and Z^ff f°r which we give only the rows corresponding to the 
evanescent operators Q25-Q28- Our results are 

/ 



7(2,0) 

■^25-28,P 



3908 
9 

1760 
3 





and 



V 



/ 



7(2,0) 

J 25-28,7V 



V 



2656 
27 


7202 
243 


157 
243 




722 
243 


55 
81 




1006 
243 




761 
162 




11302 
729 







3584 
9 


1616 
81 


3736 
81 




176 
81 


110 
27 




2192 
81 




454 
27 




4640 
243 










442528 
9 


589928 
27 




27680 
9 


15560 
9 




9344 
9 




18524 
3 




23456 
81 










145528 


860864 




41276 


108475 




2416 




177 




152512 







9 


81 




27 


54 




3 




2 




243 




11392 


3101 


4 


4 


376 


? 


? 


? 


? 


? ? 


? 


? 


? 


\ 


729 


486 


9 


" 27 


243 




















4640 
243 


949 
81 


8 
3 


8 



752 
81 


? 


? 


? 


? 


? ? 


? 


? 


? 




23456 
81 


35120 
27 


3712 

3 


04 
" 9 


1088 
9 


? 


? 


? 


? 


? ? 


? 


? 


? 




152512 
243 


337192 
81 


56 
9 


1768 
27 


688 
9 


? 


? 


? 


? 


? ? 


? 


? 


? 


/ 



(45) 



(46) 



where the question marks correspond to entries that we have not calculated. Notice that 
only the 4x4 block to the right of -^25-28 p ls nee ded to determine the O(a^) mixing of 
Qi~Qg into Qy-Qio- On the other hand the 6x4 block to the left is necessary to find the 
three-loop self-mixing of Qi-Q§ which we shall present elsewhere [23- Finally the mixing 
of evanescent operators into EOM- vanishing ones, described by Z^) 2 % N , is not necessary, 
but given for completeness here. 

Since Z^ 2 ' 2 ^ is completely determined by the one-loop mixing, we give only the non- 
vanishing building blocks of Z^ 2 ' l \ namely 

/ MI _515 _ 353 1567 67 35 58 167 64 f) \ 

36 54 486 648 243 648 w 









K2,l) 

j pp 



317 
36 


515 

54 


353 
243 


1567 
972 


07 

480 


35 
648 


58 
243 


167 

048 


64 
729 


340 
12 


3 


104 
81 


338 
81 


14 

81 


35 
108 


116 
81 


10 
27 


776 
243 








1117 
81 


31469 
324 


100 
81 


3373 
432 


10 
81 


92 
27 


1688 
243 








4079 
486 


59399 
972 


269 
1944 


12899 
2592 


50 
243 


1409 
648 


548 
729 








83080 


159926 


8839 


14573 


464 


3407 


31376 


81 


81 


81 


54 


81 


27 


243 








70100 

243 


231956 
243 


11501 

243 


78089 
648 


836 
243 


1081 
81 


21128 
729 




















650 
27 


























548 
81 


1975 
108 





























58 
3 




































58 
3 



(47) 



K2,l) 
J PN 



( 



64 
729 


3671 
1944 


1 




1 

27 


22 
243 


65 
2592 


1 

96 


1 

32 


1 

72 


1 

192 


1 

72 





1 

48 


1 

04 


776 
243 


1889 
324 


2 

3 


2 




44 
81 


259 
432 


1 

Hi 


3 
16 


1 

12 


1 

32 


1 

12 





1 

8 


3 

32 


1688 
243 


1351 
162 


16 
3 


4 



20 
81 


259 
216 


-1 


3 
8 


1 



1 

10 


1 
6 





1 

4 


3 
16 


548 
729 


3559 
972 


22 
9 


32 
27 


130 
243 


625 
648 


67 
48 


11 
4 


26 
9 


11 
24 


26 






10 
6 


11 

8 


31376 
243 


2054 
81 


304 
3 


32 
3 


320 
81 


950 
27 


35 
2 


-3 


4 

3 


1 

2 


4 
3 





-10 


3 
2 


21128 


32669 


148 


4 


1684 


12817 


35 


131 


259 


131 


250 





125 


131 


729 


486 


9 


9 


243 


648 


24 


8 


18 


48 


18 


12 


16 


















































995 7 1 
18 9 s 











43 27 c 
8 32 s 












































































































(48) 



21 



and 



7(2,1) 

J PE 



4493 

864 


49 
648 








1 

384 


35 
864 










1031 
144 


8 
9 








35 
192 
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Similarly to what happens in the case of Zp^ not all entries of Zp^ 



(49) 



are needed to find 

the 0{al) ADM of physical operators considered in this article. Needless to say, the 
mixing of physical into EOM- vanishing operators, described by Z p ^ , is not required to 
determine the mixing of physical operators at the three-loop level. However, some entries 
are important to verify the 0(a 2 s ) mixing of magnetic into non-physical operators which 
has been discussed in part in [TTj . 

As far as the mixing among EOM-vanishing operators is concerned, we have calculated 

~ (2 1) 

only the first two rows of the corresponding matrix Z N ' N . We find 
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where the question marks stand for entries that we have not calculated. 

In the case of the mixing of evanescent into other operators, we have calculated only 
the first four rows of the corresponding matrices Z^p , Z%^ an d Z^\ We get 
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and 
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Here once again question marks denote entries that we have not computed. Clearly, the 
mixing of evanescent into other operators does not affect the 0(a 3 ) mixing of physical 
operators at all, and thus is given here only for completeness. 

At the three-loop level we have calculated only a small subset of entries of Z^ 3 ' 1 ^ which 
are summarized below. Again Z^ and Z^ can in principle be obtained using Eqs. ()25j) . 
The single poles we have calculated read 
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and 
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